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I. INTRODUCTION 

In four-dimensional spacetimes like the Kerr metric [l[ , 
the existence of conserved quantities for geodesies (con- 
stants of motion) Q and the tensorial structures that 
generate them (Killing vectors, Killing tensors and 
Killing- Yano tensors ]!, H, 0] ) have been very important, 
not only elucidating particle motion in these spacetimes, 
but also leading to the separation of the Klein-Gordon 
, massless neutrino [1, 0] , massive Dira c fToL [lil ] , elec- 
tromagnetic [lH, and gravitational wave jl2|| equations. 

With the recent interest in higher-dimensional space- 
times, it has become of interest to extend these old four- 
dimensional studies to higher dimensions D. For exam- 
ple, it has been found that the rotating black hole met- 
rics [U, EH EH EH E3) El] m higher dimensions have a 
Killing- Yano tensor of rank D — 2 19, 20] (which we shall 
call a principal Killing- Yano tensor) that was used (along 
with the Killing vectors) to show [2l|, [22J that geodesic 
motion in the general D-dimensional Kerr-NUT-AdS ro- 
tating black hole spacetime [l8| is completely integrable, 
with D independent constants of motion in involution. 

For convenience, we use square brackets to denote the 
integer part of what is inside and define n ee [D/2], 
k ee [{D + l)/2], and s ee k — n (0 for even D and 1 
for odd D), so D = 2n + e = 2k - e = k + n. Then 
the Kerr-NUT-AdS spacetimes have k Killing vectors 
(giving constants of motion linear in the velocity) and n 
independent Killing tensors of higher rank (including the 
metric) [2l|, [22[ (giving other independent constants of 
motion in involution that are higher-order polynomials 
in the velocity). 

Here we show how to construct k Killing- Yano tensors, 



of ranks D — 2 j for < j < k — 1, for any spacetime 
with a principal Killing- Yano tensor. Contractions of 
each of these with itself (leaving two indices free) give k 
rank-2 Killing tensors and hence k independent constants 
of geodesic motion in involution for any spacetime with 
a principal Killing- Yano tensor. 

For the case of the Kerr-NUT-AdS spacetimes [l8[ , all 
of the k Killing vectors can also be constructed from the 
principal Killing- Yano tensor and its covariant deriva- 
tive, so all D constants of motion arise from one single 
Killing- Yano tensor (and the metric, of course, for defin- 
ing covariant derivatives and contractions). For these 
metrics, we exhibit explicitly the resulting k Killing Yano 
tensors of rank D — 2j, the k rank-2 Killing tensors, 
and the k Killing vectors. (For odd D, e = 1, one of 
these rank-2 Killing tensors is the tensor product of a 
Killing vector with itself and so is not independent or 
irreducible, leaving only D = 2k — e independent rank-2 
and rank-1 Killing tensors) We also show the relations 
of the constants of motion arising from all these Killing 
tensors with those given in [U [22j , as well as with the 
constants of motion arising from the recent separation of 
the Hamilton-Jacobi and Klein-Gordon equations [23| . 

II. GENERATING FUNCTION 

The following construction of the Killing tensors, 
Killing- Yano tensors, and of the corresponding conserved 
quantities applies for any metric [24| with a rank-2 
closed conformal Killing- Yano tensor h (or, equivalently, 
a (D — 2)-rank Killing- Yano tensor f = * h; see Sec- 
tion lIVp . Since h and / play an important role in our 
construction, we call both of them principal tensors. 

Let us recall that a rank-2 conformal Killing- Yano ten- 
sor h is an antisymmetric 2-form which obeys 

V( a /l b)c = - (gafcV e fe C c - V e /i e ( a 5fc) c ) . (1) 



'Electronic address: Pavel.Krtous@mff.cuni.cz 

^Electronic address: kubiznak@phys.ualberta.ca 

^Electronic address: don@phys.ualberta.ca 

§ Electronic address: frolov@phys.ualberta.ca 



Assuming the existence of such a tensor, we define the 



2 



2-form F = u ■ * (u ■ * h) , which in components reads 
Fab = whab- u a u c h c b - h ac U C Ub = w h c d P c a P" b ■ (2) 

Here u a is the velocity [25(, w — u a u a , and 
P c a = S c a — w~ 1 u c u a is the projector to the space 
orthogonal to u. The 2-form F is covariantly conserved 
in the direction of u, 

u a V a Fbc = . (3) 
Now we introduce the generating function W(0), 

W{0) = det(7 + y/pw^F) , (4) 

where we take F (and similarly h, P, p, and Q below) to 
be the matrix of components F a b of the 2-form F. Due to 
the antisymmetry of F and properties of the determinant, 
W{0) can be rewritten as a function of j3 instead of y/]3, 
and in terms of h and P instead of F , 

W(J3) = det 1 /^/ - pw- 2 F 2 ) = det(J - \/fihP) . (5) 

Because it is constructed only in terms of covariantly 
conserved quantities F and w, the generating function is 
conserved along geodesies, and the same is true for its 
derivatives with respect to f3. We can thus define con- 
stants of motion Cj as the coefficients in the /3-cxpansion 
ofW(/3): 



OO 



(6) 



It turns out that all terms with j > n are zero. 

To evaluate the observables Cj, we first split W{0) in 
the following way: 

W{0) = W Q (f3) E(/3) , (7) 

with 

W (p) = det(l - y/0h) , 

EC9) = det(/ +7 ^p) =tr((J-^r 1 p) • ^ 

Here p a b — w~ 1 u a Ub is the projector into the u direction, 
and we used the fact that the matrix in the determinant 
in the expression for £(/3) differs from / only in the one- 
dimensional subspace given by it. The generating func- 
tion thus splits into a part Wq(/3) independent of u and 
a part £(/3) linear in p. Using the antisymmetry of h, we 
can rewrite Wo(/3) and E(/3) in terms of the conformal 
Killing tensor with components Q a b = —h a c h c b, 



W (/3) = det 1 / 2 (j + /3Q) , 

OO 

£(/?) = tr((I + PQ)~ l p) = tr(Qip) ? 

3=0 



(9) 



We shall assume that h is non-degenerate with differ- 
ent eigenvalues. This means that there exist n uniquely 



defined 2-dimensional subspaces labeled by the index 
fi = 1, . . . , n, each of which can be spanned by a pair 
of the orthonormal vectors e M and (/x = n + /i), 
in odd number of dimensions complemented with a 
one-dimensional subspace spanned by the vector 
(0 = 2n + 1), with non-zero different eigenvalues such 
that h has the form 



(10) 



Here, = e M A are mutually orthogonal 2-forms as- 
sociated with the 2-dimensional planes. 

Using this frame we can write also the conformal 
Killing tensor Q as 



<2 = E 



(11) 



Now we can write down the functions Wq and E in 
terms of the eigenvalues . For the part independent of 
u, we get 

n n 

wm = n (i + K) = E A(j) v j > ( 12 ) 

where (cf. [3, though here in a more general situation) 
£ (13) 



vi<---<i'j 



Similarly, 



/i=i 



_. oo n 

^B+E(-wek+^ 



(14) 



.23 



j=0 



with Up, up,, and being components of u with respect 

to the dual frame e^e^e . Recall that e = in even 
dimensions D = 2n and e = 1 in odd dimensions D = 
2ti+1. 

The original generating function reads 

, n j 

= ^ElE^ 1 )^' -0 ^)^ ' (15) 

where 

wi=wtr(Q l p) = uQ l ■u = u ai Q ai a2 ...Q a ' a u a (16) 

are quantities quadratic in the velocity u given by the 
l-th. matrix power Q 1 of the conformal Killing tensor Q. 
Clearly, wq — w, and Wj — J2^( u fi + u 'p,) x 'fi f° r 3 > 0- I n 
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terms of the eigenvalues x^, from the product of Eqs. ([12 
and (1141) we obtain 



w w = ^E(^ + E4 3, K+4))^' , (17) 

3=0 /u=l 

where we have introduced the quantities (cf. fl8l ]) 

A$ = J2 (18) 



f 1 < • • • <Vj 



III. CONSTANTS OF MOTION AND RANK-2 
KILLING TENSORS 

Comparing Eqs. (fTS)) and (fT7)) with Eq. ©, we can 
identify the k — n + e = [(D + l)/2] conserved quantities 
Cj (constants of geodesic motion, j = 0, . . . , k — 1), 



(=0 



(19) 

These constants are quadratic in the velocities. They can 
be generated 0] by rank-2 Killing tensors as 



(20) 



where 



1=0 



(21) 



= sA^e°e Q + A<7> (e"e" + e^e*) 



The matrix power Q l of the tensor Q is defined in 
Eq. (fl6|) . The Killing tensors are completely symmet- 
ric tensors obeying the equations 



v (o irgj = o 



(22) 



Let us remark that the constant c„ present in an odd 
number of spacetime dimensions is the square of the con- 
stant corresponding to the Killing vector cf. Sec- 
tions HV| and El 

The relation (fl9|) can be inverted using the identities 
[6)) and (|A7|) from the Appendix. We obtain 



(23) 



3=0 



and, in an odd number of dimensions, 

2 _ c " 

Here, the quantity is defined as (cf. [TJI]) 



(24) 



(25) 



The coefficients A® of the /3 -expansion of Wo(/3) are 
the sums of all different products of j different eigenvalues 
of h, cf. Eq. (fl"3"|) . Such combinations can be obtained 
taking first the j-th wedge-power of the 2-form h, 



h^= h A • • • A h = jl 



!/!<•••<!/,• 



c, yj ui Vl A • • • A U v i , 
(26) 



and contracting it with itself in all tensor indices. Indeed, 



1 



(2^j!) 5 



(27) 



where • denotes the complete contraction, i.e., B • _B = 
B a bc...B abc " , and where we used the orthogonality 
a;'' • a;" = 2<5 M ", along with the normalization (|A4j) . 

Observing that the relation between the constants 
w^^-Cj and the tensor w~ 1 F is the same as between A^ 
and h [cf. Eqs. (gj), ©, 0, ©, and l[T3])], we obtain a 
new simple expression for the constants of motion Cj , 



1 



, } l-2j pAj % pAj 



(28) 



3 mm 2 

If one defines a 1-form v with components 

v a = h ab u b , (29) 

orthogonal to the velocity 1-form it, then Eq. (J5J implies 
that 

F = u)/i + «A». (30) 
Since uA«AmA« = 0, we have 

F AJ = i^'h Aj + juP^u AvA h A ( j -V . (31) 
The total contraction of this with itself is 

pM . pAj = ir 2, }l j . h Aj 

+ 2jw 2j - 1 h Aj »(«Ai)A ft A(j - 1) ) (32) 
+ jV^^ADAfi^" 1 ') • (ttADA^ 3 '" 1 ') . 

In the total contraction of the last term, u ■ v = 0, 
and any term with a contraction of u with any of the 
fa's gives another v which combines with the original v 
to give zero by the antisymmetry of the wedge product. 
Therefore, the only nonzero parts of the last term have 
another w = u u factor, giving a total factor of w 2 ^ 1 for 
the total contraction of that term. Upon the substitution 
of Eq. ([32]) into ([28]) . the dependence on w cancels out, 
and we recover the quadratic dependence on u which 
enters through u and v. Comparing with (|20p . we can 
write the tensorial relation between K^> and h, which 
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in components reads 



(2^1)2 \° b 11 




- 4j /i [Qbl . . . 


/i" A ' l V 1 -''«A-] 


+ 2jh a[bl ... 


(33) 






-2j/i[ Qbl ... 


h a 6 /'/,, .../'„ 6 ) , 



where we have employed the definition , the identities 
(|A2[) and (|A3D . and the normalization (|A4[) . 

Recently |2l| there have been found different conserved 
quantities, 



C, 



w 



tr[(- F y] , 



(34) 



which are, however, not quadratic in velocities. Now we 
show that also these observables can be generated from 
the generating function W((3). Taking the logarithm of 
Eq. ^ and expanding it into a power series, we obtain 



\ogW(f3) = ^tr\og(l~(3w- 2 F 2 ) 



3=0 



V 



w'J 



The constants Cj are thus (up to constant factors and 
powers of w) given by derivatives of log W(/3). We also 
obtained the relation between both sets of constants 
which can be formulated as 



^ w l Cjfi> = expO — — rC, 



3=0 



j=0 



2j vji 



(36) 



Comparing different orders of (3 we get for the first four 
constants 



Cl = 


f" 




wc 2 — 


-~c 2 - 




w 2 c 3 = 




o 


W 3 C4 — 


-ICH 





(37) 



1 

— ( 

48 



1 

- — ( 

32 



1 ( 

384' 



It is shown in [22| that the observables Cj Poisson com- 
mute between each other. The relation (|36[) . which shows 
that the Cj 's are polynomial combinations of the Cj 's and 
w with constant coefficients, thus proves that also the ob- 



servables Cj are in involution, 



{ci,c 3 } = 



(38) 



This gives non-trivial relations for the corresponding 
Killing tensors, namely 



be) 



K (l ) V e K\ 3 \ = 

e(a be) 



(39) 



IV. KILLING- YANO TENSORS 

The existence of the closed conformal Killing- Yano 
tensor h guarantees the existence of the Killing- 
Yano tensor / which is obtained by the Hodge 
dual / = *h 26] . This principal Killing- Yano ten- 
sor enables one to construct a rank- 2 Killing tensor 
K a b = f aei ' eD ~ 3 fbe 1 ...e D - 3 - Here we demonstrate that 
all the Killing tensors found in the previous section can 
be constructed in a similar way. 

First, let us recall that a conformal Killing- Yano tensor 
(CKYT) [U H, HH of a general rank r is an antisym- 
metric r-form / the covariant derivative of which can be 
split into the antisymmetric and divergence parts 



v/ = ^v/ + rv/. 



(40) 



Here A is the standard anti-symmetrization and T is the 
projection onto the 'trace' part of the tensor of rank r + 1 
which is antisymmetric in the last r indices, 



t 

T A aai a — — ; T 3a\a±A \e\a2...a r ] 

D — r + 1 



(41) 



The operation T satisfies T 2 = T and TA = AT = 0. 
This means that a tensor A satisfies TA = A if and only 
if it has the form A aai ... ar = ff [ 01 a oa ... 0r ]. The diver- 
gence part T V/ thus depends only on the divergence 
V e / e ab...- The condition (|4*U|) implies that V/ does not 
have a harmonic part [3(| (given by the complement of 
the A and T projectors), i.e., / does not have a part for 
which both d f and V-/ vanishes. 

A CKYT transforms into a CKYT under Hodge du- 
ality. The antisymmetric part A\7 f transforms into the 
divergence part T V*/ and vice versa. 

A Killing- Yano tensor / is such a CKYT for which the 
divergence part is missing, i.e., 



(42) 



The dual of a Killing- Yano tensor is a closed CKYT (see 
also [26]), i.e., an r-form obeying 



v/ = rv/ 



(43) 



The wedge product of two closed CKYTs is again a 
closed CKYT [see Eq. (|A"5]l in the Appendix]. We can 
thus start with the principal closed CKYT h and con- 
struct its wedge powers ft. A - J (J = 0, . . . , k — 1), which are 
again closed CKYTs. Their duals, 



/Cj) = * , 



(44) 



are then Killing- Yano tensors of rank D — 2j. Their com- 
ponents are 

fti) —or?? e i- e 2ifj h (ak) 

Jai.-a.D-2j — t ai...a_D_ 2 j ll e 1 e 2 ■ ■ ■ ' L e2j-ie 2 j i V*") 

where e ai ...a D are components of the Levi-Civita tensor 
e and the normalization (|A4|) has been employed. 
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Now we show that these Killing- Yano tensors gener- 
ate the rank- 2 Killing tensors K^' constructed above. 
Namely, using the identity (|Alj) . we write 



(Z>-2j-l)!C?!) a 



f(j) aei...e D _ 2 j f(j) _ 
J Jbei...er>-2j 



+ S [a h a lbl 



(46) 



.h a J 



L aib 2 ■ ■ ■ "'a j b j ] 



With the help of Eq. (IA2p , we see that the last expression 
coincides with the formula (|3"3")1 for the Killing tensors. So 
we have 



(D-2j-l)!(j!)2 



i-(j) ae 1 ...e D - 2 j Aj) 
J •> be±...eD-2j 



(47) 

Starting from the principal closed CKYT h, we 
build the sequence of closed CKYTs h Aj , given explic- 
itly in Eq. (|26p . which generates Killing- Yano tensors 
= * h Aj '. These Killing- Yano tensors can be used to 
construct the rank- 2 Killing tensors given by the 

formula ijiTjl. or explicitly by Eq. or ([55]). 

In particular, in an odd number of spacetime dimen- 
sions, the last Killing- Yano tensor /( n ) cx V eg is 
a Killing vector. Obviously, the corresponding Killing 
tensor K^ n > oc is reducible. 



V. KERR-NUT-ADS SPACETIMES 

We shall now demonstrate that the structure explored 
above is fully realized in the Kerr-NUT-AdS spacctimes 
0. 

In the notation of previous sections (i.e., using the base 
of 1-forms e a = {e^, e M , e }, fi = l ) ...,n, /i = |U + n, 
= 2ra+ 1, with the 1-form e° present only for odd D), 
the Kerr-NUT-AdS metric may be written in the or- 
thonormal form 

D n 

where the orthonormal basis one-forms are 



3=0 



(49) 



The quantities A^), in terms of coordinates x^ 
are of the form of Eqs. (JT3J) and (JTHJ), <3 M = X^/U^ with 
given by Eq. ([55]), c = fltl a? , and 



^ = (-1) 1 



1 + Xx 

rr2e 



2 fe-1 



rj( s 2 -4) + 2M iU (-^) 



1-6 



3=1 



(50) 



The constants (M^,^) are related to the mass, NUT 
parameters, and angular momenta, and A is proportional 
to the cosmological constant 18]. The metric represents 
an Einstein space obeying the Einstein equation 



R ab = (D- l)Xg ab 



(51) 



Using the identities (|A6|) and (|A7j) . we find the dual 
vectors 



e - 1/2 d 



e 6 =(- C AW)- 1 % B) 
and the corresponding inverse relations 



d - Q~ 1/2 e 



A* =^' ) (-^)) 1/2e S + E W^A, (53) 



av,„ = (-MW) 1/2 e 6 . 



It is possible to prove [3l| that that the ft. found in 
[20I ] (there called fc) in this metric in all dimensions D is 
a principal closed CKYT, which in the frame (j49|) takes 
exactly the form (flTJ|) . This means that the generally 
defined eigenvalues x^ of the principal CKYT h (fT0|) co- 
incide with the chosen ('natural') coordinates x^ of the 
Kerr-NUT-AdS metric. 

We now demonstrate that from the very existence of 
this tensor one can extract all the constants of geodesic 
motion for the Kerr-NUT-AdS metrics. Namely, besides 
the k constants of motion (fT9|) connected with the rank-2 
Killing tensors (|2Tj) , also all the k isometries follow from 
the existence of the principal CKYT h. 

First of all, it was proved in [32| that in an Einstein 
space, obeying Eq. (|5Tj) . the divergence £ of a CKYT h, 



e = 



1 



D - 1 



-Vbh" 



(54) 



is a Killing vector. In particular we find £ = d^ |20| . 

Next, using Eqs. (|2"Tj) and fix?)) , we can recover n — 1 
other Killing vectors , 



(d^r=K^%e, j=i,. 



1 



(55) 



For a similar construction in 4D see [33J, [3J, [35J . Finally, 
as mentioned in Section IIV1 in odd dimensions the last 
Killing vector is given by the n-th Killing- Yano tensor 
which in the present example turns out to be d^ n . 
It would be very interesting to find under what general 
conditions on the CKYT h, and possibly on the curva- 
ture, this construction gives all the isometries present in 
the spacetime. 
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The constants of geodesic motion in the higher- 
dimensional Kerr-NUT-AdS spacetime are now com- 
pletely determined. Denoting the constants from the 
Killing vectors as 

b 3 = {d^) a u a , j = 0,...,fc-l, (56) 

we first find the frame components of velocity u^, and 
possibly uq. It follows from Eq. ([52]) that 

u, = - i J2(-xir-^ bj , u 6 = (-cA^y 1/2 b n . 

(57) 

Comparing with (|24|) we find 




which clearly illustrates the fact that the constant c„ cor- 
responds to the reducible Killing tensor K (") oc f( n )f( n ). 
The remaining components of velocity, u fll are given (up 
to signs) in terms of the constants bj, j = 0, . . . , k— 1, and 
Cj, j = 0, . . . , n — 1 [which correspond to the irreducible 
Killing tensors (2IJ|] by Eq. ([215]) . 

The existence of n rank- 2 irreducible Killing tensors 
K^> and k = D — n Killing vectors d^ Jj is closely related 
to the question of separability of the Hamilton-Jacobi 
and Klein-Gordon equations. It is shown in [3(| [37| 
that if the spacetime possesses such tensors that satisfy 
the condition C^K® = with {8.^,8^} = 0— 
which holds in our case — then there exists a so-called 
separability structure. This structure guarantees the sep- 
arability of the Hamilton-Jacobi equation and for Ein- 
stein spaces also the separability of the Klein-Gordon 
equation. 

The separability of these equations was explicitly 
demonstrated recently 23]. It turns out that the in- 
tegration constants obtained by the separation of the 
Hamilton-Jacobi equation are the quantities Cj given by 
Eq. (jT9|) . Indeed, if we transform the tetrad components 
m m , up, and Uq into the coordinate frame we find that the 
expression U^u 2 ^ + u?) corresponds to the quantity F M 
of [H[. Comparing Eqs. ([23]) and (14) of [23|, we can 
identify the constants defined above with those from [23| 
(cf. also Eqs. ([24]) and ([58]) with Eq. (15) of f^). 

VI. DISCUSSION 

We have seen that the existence of a principal Killing- 
Yano tensor, one of rank D — 2, guarantees the existence 
of k — D — n Killing- Yano tensors of rank (D — 2j), 
j = 0, . . . , k — 1, and that each of these Killing- Yano 
tensors generates a Killing tensor of rank 2, n of which 
are irreducible. In the case of the Kerr-NUT-AdS space- 
times, the principal Killing- Yano tensor also generates all 
the k Killing vectors, and hence all D of the independent 
constants of geodesic motion. 



Our results raise various questions. For example, is the 
construction of Killing vectors by a rank-(Z) — 2) Killing 
Yano tensor general, or specific to only certain metrics? 
For what classes of spacetimes are there rank-(Z? — 2) 
Killing- Yano tensors? For what classes of such space- 
times does the rank-(D — 2) Killing- Yano tensor generate 
enough Killing tensors to give D independent constants 
of motion? Are there any new Einstein metrics within 
these classes? Do these structures enable one to separate 
the Dirac, electromagnetic, and gravitational wave equa- 
tions in the Kerr-NUT-AdS spacetimes and/or in any 
other possible members of these classes? What is the re- 
lation of the existence of principal Killing- Yano tensors 
to the algebraic type of the metric? We suspect that our 
observations may be the tip of an iceberg of important 
new relations for higher dimensional spacetime metrics. 
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APPENDIX A: USEFUL IDENTITIES 

In this Appendix we list some identities used in the 
main text. Their proofs are mostly straightforward but 
also lengthy and cumbersome. 

First we list three identities for antisymmetric tensors. 
The Levi-Civita tensor satisfies 

e ai - a ^ +1 - CD e bl ... brCr+1 ... CD = r\(D - r)!^ 1 . . . ^ , 

(Al) 

the projector on the antisymmetric tensors can be split 
as 

(r + 1) f gag . . . szi = ■ ■ ■ - r s k <i ■ ■ ■ C] 

(A2) 

and finally 

h[ ai bi ■ ■ ■ h aj b 3 ] = h ai [ bl . . . h a . bj ] , (A3) 

which holds for any antisymmetric tensor h ab - 

In our convention the wedge product is proportional to 
the anti-symmetrization of the tensor product which for 
the J-th wedge-power of a 2-form h gives 

hAj = mi Ahj - (A4) 
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Next we want to show that the wedge product of 
two closed CKYTs p and q of rank r and s respec- 
tively is again a closed CKYT. It is clear that p A q is 
closed. Rewriting the wedge product p A q as the anti- 
symmetrization of their tensor product, we get 

V e (P[ 6...9cd...]) = (^eP[ab...) Qcd...] + P[ab... (V| e | 5cd... ] ) 
= 5e[a(s^+l( V lffl^6-)9 , cd...] (A5) 

+ (-l) r c4+I^...c(V| g |g^... ] )), 

where we used the property (|43|) of p and q. We can see 
by inspection that the result has the form g a [a 1 a a 2 ...a r ] > 
so, as we discussed after Eq. (l41"j) . T(p A g) = p A q, and 
hence p A q is a closed CKYT. 



If we understand Afi' with ^ = 1, ...,n and 
j = 0, . . . , n — 1 as an n x n matrix, its inverse is 
Bfa = (-xl) n - 1 - j /Uf, with Up defined in Eq. This 
means 

j=0 ^ /x=l p 

(A6) 

We also have its 'extension' for j = n, I = 0, . . . , n — 1: 
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